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Abstract 

The discrete Fourier transform of the greatest common divisor 

m 

id[o](m) = ^gcd(fc, m)a„ , 

k=l 

with a m a primitive m-th root of unity, is a multiplicative function that 
generalises both the gcd-sum function and Euler's totient function. On the 
one hand it is the Dirichlet convolution of the identity with Ramanujan's 
sum, id [a] = id * c.(a), and on the other hand it can be written as a 
generalised convolution product, id [a] = id * a <f). 

We show that id[a](m) counts the number of elements in the set of 
ordered pairs (i, j) such that i-j = a mod m. Furthermore we generalise a 
dozen known identities for the totient function, to identities which involve 
the discrete Fourier transform of the greatest common divisor, including 
its partial sums, and its Lambert series. 



1 Introduction 

In 115] discrete Fourier transforms of functions of the greatest common divisor 
were studied, i.e. 

m 

h[a}(m) = ^2h(gcd(k,m))a k ^, 

k=l 

where a m is a primitive m-th root of unity. The main result in that paper is 
the identity^] h [a] = h * c.(a), where * denoted Dirichlet convolution, i.e. 

h[a][m) =^2h(—)c d (a), (1) 

d\m 



1 Similar results in the context of r-even function were obtained earlier, see 1101 for details. 



and 

m 

c m (a)= £ a£ (2) 
fe=i 

gcd(fc,m) = l 

denotes Ramanujan's sum. Ramanujan's sum generalises both Euler's totient 
function cj) — c.(0) and the Mobius function fj, = c.(l). Thus, identity ([I]) 
generalizes the formula 

m 

5>(gcd(k,m)) = (fc*Mm), (3) 
fe=i 

already known to Cesaro in 1885. The formula shows that /i[a] is mul- 
tiplicative if h is multiplicative (because c. (a) is multiplicative and Dirichlet 
convolution preserves multiplicativity) . 

Here we will take h = id to be the identity function (of the greatest com- 
mon divisor) and study its Fourier transform. Obviously, as id : id(n) = n is 
multiplicative, the function id [a] is multiplicative, for all a. Two special cases 
are well-known. Taking a = 0we have id[0] = V, where 

in 

■p{m) = ^gcd(fc,m). (4) 

k=l 

is known as Pillai's arithmetical function or the gcd-sum function. Secondly, by 
taking a = 1 in we find that id[l] = id * /i equals <p, by Mobius inversion 
of Euler's identity <fi * u — id, where u = is the unit function defined by 
u(m) = 1. 

Let T™ denote the set of ordered pairs (i, j) such that i ■ j = a mod in, the 
set of factorizations of a modulo m. We claim that id[et](m) counts its number 
of elements. Let us consider the mentioned special cases. 

a = For given i G {1,2,..., to} the congruence i • j = mod to yields 

— — rj = mod — — r, 

gcd(z,mj gcd(z,TO) 

which has a unique solution modulo m/gcd(j, to), and so there are gcd(i, to) 
solutions modulo to. Hence, the total number of elements in J 7 ™ is V(ra). 

a = 1 The totient function <j>(m) counts the number of invertible congruence 
classes modulo to. As for every invertible congruence class i modulo to 
there is a unique j = i^ 1 mod to such that i ■ j = 1 mod to, it counts 
the number of elements in the set J 7 ™. 

To prove the general case we employ a Kluyver-like formula for id [a], that is, a 
formula similar to the formula for the Ramanujan sum function 

d\gcd(a.k) 
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attributed to Kluyver. Together the identities ([I]) and ^ imply, cf. section |3j 

Cd[a](m)= Yl ( 6 ) 

<i|gcd(a,m) 

and we will show, in the next section, that the number of factorizations of a 
mod m is given by the same sum. 

The right hand sides of ([5| and ^ are particular instances of the so called 
generalized Ramanujan sums pQ, and both formulas follow as consequence of a 
general formula for the Fourier coefficients of these generalised Ramanujan sums 
[H [3] . In section [3] we provide simple proofs for some of the nice properties of 
these sums. In particular we interpret the sums as a generalization of Dirichlet 
convolution. This interpretation lies at the heart of many of the generalised 
totient identities we establish in section [H 



2 The number of factorizations of a mod m 

For given i, m G N, denote g = gcd(i,m). If the congruence i ■ j = a mod m 
has a solution j, then g \ a and j = i~ l a/g is unique mod m/ g, so mod m there 
are g solutions. This yields 

m 

»=1 

god(4,m)|o 



which can be written as 

m 

#-^ n = E E d ( ? ) 

d\a i=l 

gcd(i,m) = d 

If d \ m then the sum 

m 

E i 

i=l 

S cd(i, m )=d 

is empty. Now let d \ m. The only integers i which contribute to the sum are 
the multiples of d, kd, where gcd(k,m/d) = 1. There are exactly 4>{m/d) of 
them. Therefore the right hand sides of formulae ^ and ([7| agree, and hence 
#77 = id[a](m). 



3 A historical remark, and generalised Ramanu- 
jan sums 

It is well known that Ramanujan was not the first who considered the sum 
c m (a). Kluyver proved his formula ^ in 1906, twelve years before Ramanujan 
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published the novel idea of expressing arithmetical functions in the form of a 
series ^ s a s c s (n) [IT]. It is not well known that Kluyver actually showed that 
c m (a) equals Von Sterneck's function, introduced in [13], i.e. 

u( ™ , U(m) 

r v gca(a.m) ' 

This relation is referred to in the literature as Holder's relation, cf. the remark 
on page 213 in [JJ. However, Holder published this relation thirty years after 
Kluyver 7 . We refer to [IJ Theorem 2], or [31 Theorem 8.8] for a generalisation 
of pi). The so called generalized Ramanujan sums, 

f* a g(m)= Msfy, ( 9 ) 

d\gcd(a,m) 

were introduced in [JJ. The notation * a is new, the sums are denoted S(a; m) in 
PQ, s m (a) in [2], and S/^Ojm) in [3]. In the context of r-even functions [10] the 
sums are denoted S/ jff (a), and considered as sequences of m-even functions, with 
argument a. We consider the sums as a sequence of functions with argument 
m, labeled by a. We call / * a g the a-convolution of f and g. 

The concept of a-convolution is a generalization of Dirichlet convolution as 
/ *o 9 = / * 9- As we will see below, the function / * a g is multiplicative (for 
all a) if / and g are, and the following inter-associative property holds, cf. [3, 
Theorem 4]. 

if* a g)*h = f* a (g*h). (10) 

We also adopt the notation f a = id * a /, and call this the Kluyver, or a- 
extension of f . Thus, we have /o = id * /, fi = /, and formulas ^ and (|6| 
become c m (a) = /x a (m), and id [a] = <j) a , respectively. 

The identity function /, defined by /*/ = /, is given by I(k) = [k = 1], 
where the Iverson bracket is, with P a logical statement, 

\P) = l 1 ifP ' 
L J \ if not P. 

Let us consider the function / * a I. It is 

f* a I(k)= J2 f(d)[d = k] = [k\a]f(k). 

d|gcd(a,fc) 

Since the function k — > [k | a] is multiplicative, the function / * a / is multiplica- 
tive if / is multiplicative. Also, we may write, cf. [3l eq. (9)], 

/ *a g(m) = J> I a]f(d)g(^) - (/ * a 7) * g(m), 

d\m 

which shows that / * a g is multiplicative if / and g are. Also, the inter- 



associativity property ( 10 ) now easily follows from the associativity of the 
Dirichlet convolution, 

(/ * a g)*h = ((/ * a I) * g) * h = (f * a I) * (g * h) = f * Q (g * h). 
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We note that the a-convolution product is neither associative, nor commutative. 
The inter-associativity and the commutativity of Dirichlct convolution imply 
that 

f a *g = (f*g) a = f*g a - (n) 

Formula (|6| states that the Fourier transform of the greatest common divisor 
is the Kluyver extension of the totient function. We provide a simple proof. 

Proof [of ([©J)] Employing |IJ, |i| and ( |ll| we have id [a] = id*e.(a) = id * /Lt a = 
(id * n) a = <j> a . □ 

The formula ^ also follows as a special case of the following formula for the 
Fourier coefficients of a-convolutions, 



f* a g{m) = ^h k {m)a k ^, 



hk = 



• *fe 



k=l 



f_ 
id' 



(12) 



given in [HE]- The formula (12) combines with (JlJ and ^ to yield a formula 
for functions of the greatest common divisor, h[k] ■ m — > /i(gcd(fc, m)), namely 



h[k] = (h* /J,) *fc u. 



(13) 



Proof [of (13)] The Fourier coefficients of h[a](m) are h[k](m). But h[a] — 
h*c,(a) = (id* a n)*h = id* a (fj,*h), and so, using (12), the Fourier coefficients 
are also given by (h * fi) *k u(m). □ 
For a Dirichlet convolution with a Fourier transform of a function of the greatest 
common divisor we have 

f*g[a] = f*g[a]. (14) 

Proof [of (l4J)] / * g[a] = / * (g * jti ) = (/ * g) * \i a = f * g[a] □ 
Similarly, for an a-convolution with a Fourier transform of a function of the 
greatest common divisor, 



f*at[b] = f* a g[b]. 



Proof [of ( 15 1] / * a g[b] = f * a (.9 * W>) = (/ = / *o g[b] 



(15) 
□ 



4 Generalised totient identities 

The totient function is an important function in number theory, and related 
fields of mathematics. It is extensively studied, connected to many other notions 
and functions, and there exist numerous generalisation and extensions, cf. the 
chapter "The many facets of Euler's totient" in (12]. The Kluyver extension of 
the totient function is a very natural extension, and it is most surprising it has 
not been studied before. In this section we generalise a dozen known identities 
for the totient function <j), to identities which involve its Kluyver extension <f> a , 
a.k.a. the discrete Fourier transform of the greatest common divisor. This 
includes a generalisation of Euler's identity, the partial sums of <p a , and its 
Lambert series. 
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(p k -p k - l ){l + l) ifl<k, 
(k + l)p k - kp k - x if Z > fc, 



(16) 



4.1 The value of a at powers of primes 

We start by providing a formula for the value of (j> a at powers of primes. This 
depends only on the multiplicity of the prime in a. The formulae, with p prime, 

V( P k ) = (k + l)p k - kp k ~\ <p{p k )=p k -p k ~\ 

of which the first one is Theorem 2.2 in [4], generalise to 

where I is the largest integer, or infinity, such that p l \ a. 
Proof [of ([16])] We have 

Mp k ) = E ^(7) 

d\gcd(p l ,p fc ) 
min(/,fc) 

r=0 

1 (E^oV-p^+p* in > fc , 



which equals (16 1. 



4.2 Partial sums of a /id 

To generalise the totient identity 



■A ^(fc) _ A fi(k) n , 

^ fc ^ fc L fc J ' 

fe=i fe=i 



to an identity for (/) we first establish 



/o(fc) = f(k) , n 
Z^ z- Z^ t L;J 



fc=i 



(17) 



(18) 



Proof [of ( IS )] Since there are \n/d\ multiples of d in the range [l,n] it follows 
that 

/(rf).n 



fc=i 



E f_ * id(fc) = y> >p f(d) _ | 



As a corollary we obtain 



y> / *a go(fc) = y> f*a9(k) , n 



□ 



(19) 



fc=i 



fc=i 
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Proof [of ^] Employing ([To| we find 

■A f* a (g*id)(fc) = A (f* a g) *id(fc) _ A / * a g(fc) n 



fc=i 



fe=i 



k=l 



Now taking / = id and g = fi in ( 19 ) we find 

<Pa(k) 



E 



fc ^-^ fc k 

k=l k=l 



E Cfc(a) , n 



□ 



(20) 



4.3 Partial sums of V a /id expressed in terms of <p a 

Taking / = id and g = <j> in (|19| we find 



E Vajk) _ A a (fc) n , 
fc fc '-fe-'" 



fe=i 



fc=i 



(21) 



involving the totient function and the gcd-sum function, 



Note that by taking either a = in (20 1, or a = 1 in (21 1, we find an identity 

(22) 



h ~ 2^1 h L.J 



</>(fc) n 



fe=i 



fe=i 



k L k- 



4.4 Partial sums of <p a 

To generalise the totient identity, with n > 0, 



we first establish 



fe=i \ fe=i / 

n 1 / n n \ 

fe=i \fc=i fe=i / 



(23) 



(24) 



Proof [of ( 24 )] We have, by changing variable fc = dl, 



£(2/ * id - / * u)(k) = £ 1} 

fe=l fc=l d\k 

n ln/d\ 

= EE 

d=l 2 = 1 



□ 
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Note that this gives a nice proof of (23), taking / = fi, as Y2k=i ^(^) = [k > 0]- 
When / = /i a , then ( 11 ) implies / * id — <j) ai and / * u — I ai and therefore as a 



special case of ( 24 1 we obtain 



n .J / n \ 

E^( fc ) = 2 E fc t fc - n] + E c *(«)^J 2 

fe=l \fe|a fc=l / 



(25) 



We remark that when n > a we have Y)u a k[k < n] = tr(a), where a — id * u is 
the sum of divisors function. 



4.5 Generalisation of Euler's identity 

Euler's identity, <f> * u = id, generalises to 

^ (t>a{d) = r(gcd(a, m))m 

e£|m 

where r = it * u is the number of divisors function. 



Proof [of (26)] We have 4> a * u = ((/) * u) a = id a where 



id a (m) = ^ d-j = mr(gcd(a,m)). 



d|gcd(a,m) 



(26) 



(27) 
□ 



4.6 Partial sums of V a expressed in terms of <p a (and r) 



If / = <f> a , then / * id = V a , and (|24j) becomes, using (26 1 



n / n n \ 

E = 2 E r(gcd(a ' k))k + E ^( fc ) l^j 2 ■ (28) 

fc=l \fc=l k=l ) 

4.7 Four identities of Cesaro 

According to Dickson [3] the following three identities were obtained by Cesaro: 

(29) 



E#y =nn), 

d\n 

E^)^(?) = E^fe cd ^»)- 



(30) 
(31) 
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Identity (291, which is Theorem 2.3 in [3], is obtained by taking a = in ([6]), or 
h = id in (|3 1 . It generalises to 



d\n 



Proof [of ((32])] By taking / = (f> and g = id in ( flT] ). 

Identity ( 30 ) is obtained by taking /i = 1 /id in ((3] and generalises to 



d|n J=l ci|gcd(a,n) & w ' d y 



Proof [of ((33])] By taking / = <j> and g = 1/id in ([III. 

Identity (31 ) is also a special case of ((3]), with h — (f>. It generalises to 

n 

5>«o*)&(?) = E E ^(s cd (^?))- 



(32) 



□ 



(33) 



□ 



(34) 



d\n 



j = l d|gcd(a,n) 



Proof [of ([34])] We have 

(id * a (j>) * (id * b <$) = id * a (id * b (0 * 0)) = id * a (id * b <p[0]) = id * a <f> b [0], 
and evaluation at m yields 



d|gcd(a,m) j = l 



d|gcd(a,m) j=l 



The more general identity ^ generalises to 

m 

^2 /i Q (gcd(fc, m)) = h * (j} a (m). 



□ 



(35) 



k=l 



4.8 Three identities of Liouville 

Dickson [FJ p. 285-286] states, amongst many others identities that were pre- 
sented by Liouville in the series [5], the following 



5>(d)r(^)=«7(m), 



^0(d)a-[n + l](— ) = ma[n](m), 

d\m 



E^M^) = E d O 



(36) 
(37) 

(38) 
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where a[n] — id[n]*u, id[n](m) = m n , and 9(m) is the number of decompositions 
of m into two relatively prime factors. All three are of the form 4> * f = g and 
therefore they gain significance due to t hought Liouville might not have 
been aware of this. For example, ^ and (36 1 combine to yield 

m 

r(gcd(fc, m)) = a(m). 

k=i 

The three identities are easily proven by substituting r = u * u, a[n] = id[n] * u, 
t o id [2] = 8 * u, = * id, and using \i * it = /. They generalise to 



^<l>a{d)T{ — ) = cr a (m), 
^0 Q (d)cr[?i+ l](-r) = ™«*a cr[n](m), 

E Ud)r{~) = E dT(gcd(a, d))0(^). 



(39) 
(40) 

(41) 



These generalisation are proven using the same substitutions, together with 
dill, or for the latter identity, (TlOb and d27|. 



4.9 One identity of Dirichlet 

Dickson [5] writes that Dirichlet 6 , by taking partial sums on both sides of 
Eulcr's identity, obtained 



k=l 



n + l 
2 



By taking partial sums on both sides of equation ( 26 1 we obtain 



k=l 



E^w=E d 

d\a 



(42) 



Proof [of (42)] Summing the left hand side of (26) over m yields 

n n 

EE^w = EL?J^(d) 



m— 1 ( 



d=l 



and summing the right hand side of ( 26 ) over m yields 

\n/d\ 



E T(gcd( a , m))m = e E m = EE* = E d L™j (l^j + l) A 

ro= 1 d|gcd(a,m) d|a fc=l d|a 



m— 1 



□ 
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4.10 The Lambert series of a 

As shown by Liouville [9], cf. [3 p. 120], the Lambert series of the totient 
function is given by 



m— 1 

The Lambert series for cj) a is given by 



{l-xf 



T,Mrn)^^p[a]( X )—^— 2 , (43) 

m=l V ' 

where the coefficients ofp[a](a;) = 52fc=i c[a\{k)x k ~ 1 are given by c[a] = id a oi[<z], 
and t[a] is the piece-wise linear function i[a](rt) = a — \n — a\. 

The polynomials p[a] seem to be irreducible over Z and their zeros are in 
some sense close to the a-th roots of unity, see Figures [T] and [2j 



Figure 1: The roots of p[37] are depicted as crosses and the 37 th roots of unity 
as points. This figure shows that when a is prime the roots of p[a] that are close 
to 1 are closer to a th roots of unity. 

Proof [of ( 43 1] Cesaro proved [S] 

oo n oo 



n=l 



n=l 



cf. exercise 31 to chapter 2 in [Mj. By substituting (261 in this formula we find 



Z2<t>a(n) 1 _ xn =J2x n T(gcd(a,n))7 
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Multiplying the right hand side by (1 — 2x a H- x 2a ) yields 

oo oo 

(£zV(gcd(a,n))n)-2( ]T x"r(gcd(a, n)){n - a)) 

n— 1 n— a+1 

oo oo 

+ ( £ x n T(gcd(a,n))(n-2a)) = J24a](n)x n , 



where 
c[a](n) 



r(gcd(a, n))n < n < a, 

r(gcd(a, n))(n — 2(n — a)) = r(gcd(a, n))(2a — n) a < n < 2a, 
r(gcd(a, n))(n — 2(n — a) + n— 2a) =0 n > 2a. 



Rewriting, using (271, the fact that gcd(a, a + k) = gcd(a, a — fc), and dividing 
by x, yields the result. □ 



Figure 2: The roots of p[35] are depicted as crosses and the 35 th roots of unity 
as points. This figure shows that roots of p[a] are closest to primitive a th roots 
of unity. 

4.11 A series related to the Lambert series of cf) a 

Liouville [9 also showed 



V Mm)— = (1 + x 2 )^ 

We show that 



2\2 ' 



E ^'Ni^ = gMto n _ X x 2a)2 ■ ( 44 ) 

m— 1 V ■ V 
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where g[a](a;) = J2kLi b[a](k)x , with 

b[a] = h[a] o t[2a], h[a](k) = id a (fc) - 2 [2 | fc]id a (fc/2). 

At the end of this section we show that 1 + x 2 divides g[a](a;) if a is odd. 

Proof [of (44)] As the left hand side of (44 1 is obtained from the left hand side 
of ( 43 1 by subtracting twice the same series with x replaced by x 2 , the same is 
true for the right hand side. Thus it follows that g[a](a;) = p[a](x)(l + x a ) 2 — 
2p[a]{x 2 )x 1 and hence, that 



b[a](k) = 



id (fc) - 2 [2 | fc]id (fc/2) k < a, 

id a (2a - fc) + 2id a (fc - a) - 2 [2 | fc]id (fc/2) a < k < 2a, 

2id Q (3a — k) + id a (fc - 2a) - 2 [2 | fc]id a (2a - fc/2) 2a < fc < 3a, 

id a (4a - fc) - 2 [2 | fc]id Q (2a - fc/2) 3a < fc < 4a. 



The result follows due to the identities 

id a (2a - fc) + 2id a (fc - a) = id a (fc), 2id a (3a - fc) + id Q (fc - 2a) = id a (4a - fc), 



which are easily verified using (27 1. □ 

We can express the functions b[a) and h[a] in terms of an interesting fractal 
function. Let a function k of two variables be defined recursively by 



2 | n,2{a, or n = 0, 

K[a](n)={ «[a/2](n/2) 2|n,2|a, 
r(gcd(a,n)) 2\n. 



(45) 



The following properties are easily verified using the definition. We have 

K[a](2a + n) = «[o](2a - n), (46) 

and, with gcd(a, 6) = 1, 



n[an](bn) 



2 | 6, 
a(n) 2\b, 



(47) 



where a denotes the number of odd divisors function, i.e. for all fc and odd m 

a(2 k m) = r(m). (48) 



Property ( 47 1 is a quite remarkable fractal property; from the origin in every 

(49) 

b[a] = K[a]t[2a]. (50) 



direction we see either the zero sequence, or a, at different scales. 
We claim that 

h[a] = k [a] id 

follows from fl45l), (|48l, and &f\. From (|49| and (|46l we obtain 
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We now prove that 1 + x 2 divides q[a] when a is odd. Noting that b[a](2a + k) = 
b[a](2a — k) and, when 2 \ a, b[a](2k) — 0, we therefore have 

2a 

q[ a ](x) = 5^6[o](2n- l)x 2n ~ 2 

n=l 
a 

= b[a](2a -2m + l)x 2a ~ 2m + 6[o](2a + 2m - lj^Sa+am-a 

m— 1 
a 

= J] 6[a](2a - 2m + l)x 2a - 2m (l + x 4m - 2 ), 

m— 1 

which vanishes at the points where a; — — 1. □ 
Apart from the factor 1+x 2 when a is odd, the polynomial q[a] seems to be 
irreducible over Z and its zeros are in some sense close to the 2a th roots of —1 
or, to the (a + l) st roots of unity, see Figure [3] 



-0.5 



Figure 3: The roots of q[19] are depicted as boxes, the 38 th roots of —1 as 
points, and the 20 th roots of unity as crosses. 

4.12 A perfect square 

Our last identity generalises the faint fact that (f)(1) = 1. We have 

n 

£>„(n)=n 2 . (51) 



Proof [of (51 1] For any lattice point (i, j) in the square [1, n] x [1, n] the product 
i ■ j mod n is congruent to some a in the range [1, n]. □ 
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